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Abstract 

The first atmospheric Cherenkov telescope of VERITAS (the Very Energetic Ra- 
diation Imaging Telescope Array System) has been in operation since February 
2005. We present here a technical description of the instrument and a summary of 
its performance. The calibration methods are described, along with the results of 
Monte Carlo simulations of the telescope and comparisons between real and simu- 
lated data. The analysis of TeV 7-ray observations of the Crab Nebula, including 
the reconstructed energy spectrum, is shown to give results consistent with earlier 
measurements. The telescope is operating as expected and has met or exceeded all 
design specifications. 
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1 Introduction 



The first detection of Clierenkov liglit produced by cosmic ray air sliowers in 
tlie atmospliere was made over lialf a century ago [1]. Tlie detection metliod 
remains essentially unchanged today; relatively crude mirrors are used to re- 
flect Cherenkov photons onto a photomultiplier tube (PMT) detector pack- 
age, and fast electronics discriminate the brief (few nanoseconds) Cherenkov 
flashes from the background night-sky light. The flrst successful application of 
this technique to 7-ray astronomy was realised by the Whipple Collaboration, 
who used a 10 m diameter light collector and an array of PMTs to record 
images of the angular distribution of Cherenkov light produced by air showers 
[2]. The shape and orientation of the images were used to efficiently select 
7-ray initiated air shower candidate events from among the otherwise over- 
whelming background of cosmic ray initiated events [3]. The Whipple 10 m 
telescope provided the flrst detection of an astrophysical 7-ray source, the 
Crab Nebula, using this imaging technique [4]. Arrays of imaging atmospheric 
Cherenkov telescopes provide a further increase in sensitivity and in angular 
and energy resolution, as demonstrated by the HEGRA experiment [5]. The 
stereoscopic imaging atmospheric Cherenkov technique is now being exploited 
using large reflectors by four projects worldwide; MAGIC [6] and VERITAS 
[7] in the Northern hemisphere and HESS [8] and CANGAROO III [9] in the 
South. HESS was the flrst of these to come online and recent observations 
have produced a wealth of new discoveries [10]. 

The flrst stage of the VERITAS project will use four telescopes at a site in 
Arizona. Each telescope will consist of a 12 m diameter segmented reflector 
instrumented with a 499 element PMT imaging camera. The flrst of these 
telescopes has been installed at the base camp of the Whipple Observatory at 
Mt. Hopkins (at an altitude of 1275 m) and saw flrst light in February 2005. 
We begin by describing the various elements of the detector and assess their 
technical performance. The calibration and data analysis procedures are then 
discussed, and the overall performance of the instrument is characterised using 
Monte Carlo simulations and the results from preliminary observations of TeV 
7-ray sources. 



2 The Telescope 



Figure 1 shows the flrst complete VERITAS telescope installed at the Whipple 
Observatory base camp. The building closest to the telescope contains the 
trigger and data acquisition electronics and is the control room from which 
the telescope is operated. The telescope was operated initially in 2004 as a 
prototype, with one third of the mirror area and half of the PMTs of the 



completed telescope [11]. 

2.1 Mechanics and Tracking 



The basic mechanical structure of the telescope is similar to that of the Whip- 
ple 10 m telescope, consisting of an altitude-over-azimuth positioner and a 
tubular steel Optical Support Structure (OSS). The camera is supported on 
a quadropod, and a mechanical bypass of the upper quadropod arm transfers 
this load directly to the counterweight support. The positioner is a commercial 
unit manufactured by RPM-PSI (Northridge, California); the OSS is a steel 
space frame, custom designed by M3 engineering (Tucson, Arizona) and fab- 
ricated by Amber Steel (Chandler, Arizona) [12]. The maximum slew speed is 
measured to be 0.3° s~^. Tests with a modified drive system have proved that 
it is possible to reach maximum slew speeds of 1° s~^; the first telescope will 
shortly be upgraded and the remaining VERITAS telescopes will have this 
modification installed as standard. The telescope encoder measurements are 
written to a database at a rate of 4 Hz and indicate that the tracking is stable 
with a relative raw mechanical pointing accuracy of typically < ±0.01° (Fig- 
ure 2). Various methods for measuring and improving the absolute pointing 
accuracy are under investigation. 



2.2 Optics 



The telescope optics follows a Davies-Cotton design [13], but with a 12 m 
aperture reflector and a 12 m focal length. The reflector comprises 350 hexag- 
onal mirror facets (Figure 3), each with an area of 0.322 m^, giving a total 
mirror area of ~ 110 m^. The use of hexagonal facets allows the full area of 
the OSS to be exploited. The facets are made from glass, slumped and pol- 
ished by DOTI (Roundrock, Texas), then cleaned, aluminized and anodized 
at the VERITAS optical coating laboratory. The reflectivity of the anodized 
coating at normal incidence is shown as a function of wavelength in Figure 4; 
it is typically > 90 % at 320 nm. Each facet has a 24.0 m ± 1 % radius of 
curvature and is mounted on the spherical front surface of the OSS (radius 
12 m) using a triangular frame. Three adjustment screws allow each facet to 
be accurately aligned. 

The reflector facets are aligned manually using a laser system installed at a 
point facing the centre of the reflector at a distance of twice the focal length 
(24 m). An initial alignment of the mirror facets resulted in a point spread 
function (PSF) of ~ 0.09° FWHM. The 7-ray observations reported on in this 
paper were made under this alignment condition. A subsequent realignment 
reduced the PSF to 0.06° (Figure 5) at the position of Polaris (elevation 31°), 



degrading at higher elevations due to flexure of the OSS. The technique of bias 
alignment, wherein the mirror facets are aligned such that the PSF is optimum 
over the most useful observing range, has been successfully employed on the 
Whipple 10 m telescope in the past and will be applied to the VERITAS 
telescopes to achieve a PSF of < 0.06° FWHM over the 40° - 80° elevation 
range. 



2. 3 Camera 



Figure 6 shows the front face of the telescope focus box. The box is 1.8 m 
square and allows for future expansion to increase the camera field-of-view. 
The imaging camera's 499 pixels are Photonis 2.86 cm diameter, UV sensitive 
PMTs (model XP2970/02), with a quantum efficiency > 20% at 300 nm. 
The angular pixel spacing is 0.15°, giving a total field-of-view of diameter 
~ 3.5°. Light cones have not yet been installed; two different designs are 
being fabricated and will be tested. They are expected to significantly improve 
the total photon collection efficiency as well as to reduce the rate of off-axis 
background photons in each PMT. 

The PMT high voltage is provided by a multichannel modular commercial 
power supply (CAEN) which allows each PMT to be controlled individually. 
The high voltage is chosen to give a PMT gain of ~ 2 x 10^. The signals are 
amplified by a high-bandwidth pre-amplifier integrated into the PMT base 
mounts. This circuit also allows the PMT anode currents to be monitored and 
calibrated charge pulses to be injected into the signal chain. Average currents 
are typically 3 /lA (for dark fields) to 6 /iA (for bright fields), corresponding to 
a night-sky photoelectron background of 100-200 MHz per PMT at this site. 
Signals are sent via ~ 50 m of 75 f] stranded cable (RG59) to the telescope 
trigger and data acquisition electronics which is housed in the control room. 



2.4 Trigger 



In an imaging atmospheric Cherenkov telescope, precise timing between trig- 
ger channels is desirable in order to reduce the coincidence resolving time and 
hence lower the detector energy threshold. To achieve this, each channel is 
equipped with a custom-designed constant fraction discriminator (CFD) [14] 
which has a trigger time which is independent of the input pulse height. The 
CFD output width is programmable in 12 steps between 4 and 25 ns; a width 
of 10 ns was used as standard for telescope operations in 2005. A 3-bit, 6 ns 
programmable delay is provided for each channel so as to correct for system- 
atic differences in the relative signal paths due to cable length differences and 
the volt age- dependent PMT transit times. 



The CFD signals are copied and sent to a topological trigger system which is 
similar to that used successfully on the Whipple 10 m telescope [15], but with 
an improved channel-to-channel timing jitter of < 1 ns. The system contains a 
memory look-up which can be pre-programmed in a few minutes to recognise 
patterns of triggered pixels in the camera; for example, any 3 adjacent pixels. 
The required overlap time between adjacent CFD signals is ~ 6 ns. The topo- 
logical trigger system reduces the rate of triggers due to random fluctuations of 
the night-sky background light and preferentially selects compact Cherenkov 
light images. Figure 7 shows the trigger rate as a function of CFD threshold 
for two different topological trigger conflgurations. Observations in 2005 were 
all made with a rather conservative CFD threshold corresponding to ~ 6 — 7 
photoelectrons and a 3- fold adjacent pixel topological trigger conflguration, 
giving a cosmic ray trigger rate at high elevation of ~ 150 Hz. 

VERITAS will also include an array level trigger system which triggers the 
data acquisition on coincident events over a pre-deflned number of telescopes. 
A preliminary array level trigger is installed on the flrst VERITAS telescope 
which is used to latch the event number and GPS time for each event. 



2.5 Data Acquisition 



The centrepiece of the data acquisition chain is a custom-built 500 MHz flash- 
ADC system [16]. Each PMT signal is digitized with an 8-bit dynamic range 
and a memory depth of 32 /xs. By default, the signal traces follow a high gain 
path to the FADC. If the dynamic range is exceeded, an analog switch connects 
the FADC chip to a delayed low gain path instead, thus extending the dynamic 
range for each 2 ns sample. The electronic noise is small, with a sample-to- 
sample standard deviation of ~ 0.5 digital counts and an event-to-event stan- 
dard deviation over a 10 sample integration window of ~ 1.5 digital counts 
Figure 8 shows a typical FADC trace in a signal channel. 

The FADCs are deployed in four custom VME crates, where they are read out 
by local single board computers (VMIC). Buffered events from the single board 
computers are transferred via a Scaleable Coherent Interface to an event- 
building computer, where they are integrated, tested, and passed on to the 
online analysis system and data harvester. The FADC readout window size and 
position are programmable; a 24 sample readout on all 500 channels results in 
a data size of 13.5 kb per event and a deadtime of ~ 10 % at 150 Hz. While 
this is manageable for a single telescope, the VERITAS array will produce 
four times as much data at higher rates (up to 1 kHz). To cope with this, a 
zero suppression scheme has been implemented wherein only those channels 
with a peak signal larger than some preset value are read out, reducing the 
data size by a factor of > 4. The CFDs are integrated onto the FADC boards 



and the CFD hit pattern and the rate of single channel triggers can also be 
read out to the data stream. 

The FADCs give two principal benefits over simple charge integrating ADCs. 
They allow the application of digital signal processing techniques, for example, 
actively placing and minimizing the charge integration gate, thus improving 
the signal/noise per PMT and lowering the effective energy threshold in the 
analysis. In addition, they provide measurements of the time distribution of the 
Cherenkov photons across the image, which may help to reject the hadronic 
cosmic ray background showers and improve the accuracy of 7-ray shower 
parameter reconstruction [17]. 



3 Calibration 



The telescope calibration is divided into two sections. Absolute calibration 
is concerned with understanding the signal size produced by a single photon 
such that detector response and the energy scale can be accurately modelled in 
Monte Carlo simulations. Relative calibration involves measuring the various 
calibration constants for each signal channel so as to flat-field the response of 
the camera in the data analysis. 



3.1 Absolute Calibration 



The overall photon conversion factor of the telescope is a combination of the 
mirror reflectivity, the collection efficiency and quantum efficiency of the PMT 
photocathode, and the conversion factor of the electronics chain. The mirror 
reflectivity and PMT photocathode characteristics can be determined from 
laboratory measurements and manufacturer's specifications; the conversion 
factor of the electronics chain is defined as the average integral number of 
FADC counts produced by the generation of a single photoelectron in the 
PMTs and measured in situ. A nitrogen laser is used to deliver a ~ 4 ns 
wide pulse to a dye module which fluoresces at 400 nm. This light is passed 
via optical fibre to an opal glass diffuser situated in front of the camera at 
a distance of ~ 4 m, so as to provide uniform illumination [18]. The camera 
face is covered by a semi-reflective mylar sheet which shields the PMTs from 
background light, and the laser pulse intensity is adjusted using filters to 
provide an average illumination of ~ 1 photoelectron per PMT. A delayed 
copy of the laser trigger is used to trigger the data acquisition such that the 
signal arrives centered in the readout window. Figure 9 shows the histogram 
of the integrated FADC counts for a single PMT at standard operating gain. 
The conversion factor averaged over all PMTs is 0.19 ± 0.02 photoelectrons 



per digital count. 

Two complementary methods have been used to confirm the overall photon 
collection efficiency of the telescope. Single muons within the field-of-view 
produce sharply defined ring images, and the expected amount of Cherenkov 
light per unit arc length can be accurately calculated, providing a well cali- 
brated natural "test beam" [19]. A comparison of the total charge in simulated 
muons and in real muons selected from data taken with the first VERITAS 
telescope indicates that the total photon collection efficiency is well matched 
[20]. A novel alternative approach to absolute calibration, borrowed from at- 
mospheric fiuorescence detectors of ultra-high energy cosmic rays [21], has also 
been tested. A calibrated laser pulse directed towards the zenith is used to 
produce a fiash of Rayleigh-scattered light with an intensity that can be cal- 
culated very accurately when atmospheric conditions are good. Preliminary 
results with this method again indicate good agreement between simulations 
and data [22]. 



3.2 Relative Calibration 



Prior to the parameterization of Cherenkov images in the camera, each signal 
channel must be calibrated. A pulse generator is used to trigger the data 
acquisition at a rate of 3 Hz during the data-taking in order to generate events 
with no Cherenkov light present. A histogram is constructed of the integrated 
number of digital counts in each FADC trace (for a given integration window 
size); the mean of this histogram is the electronic pedestal value and the width 
is the pedestal standard deviation, which provides a measure of the night sky 
background noise level. Pedestal and pedestal standard deviation values are 
currently calculated once per 28 minute data run. The same laser system as 
used for the single photoelectron calibration is used once per night to provide 
~ 1000 bright, uniform, time-coincident photon pulses across the camera. The 
mean of a histogram of the integrated signal in each FADC trace over these 
events measures the relative gain of the signal channels. Similarly, a histogram 
of the arrival time of the laser pulse in each trace measures the relative time 
offsets for each channel. Relative PMT gains are set to within ~ 10% by 
adjusting the PMT high voltages. Approximately 95% of the signal channels 
show a relative time offset of less than ±3 ns; this is compensated for in the 
trigger hardware as discussed in Section 2.4. 



4 Data Analysis and Results 



Development work to define the optimum data analysis methods for this new 
instrument is ongoing. The data products, in particular the FADC informa- 
tion, provide the opportunity of exploiting many techniques which have not 
been applicable to previous generations of Cherenkov telescopes; for exam- 
ple, the use of digital signal processing algorithms to parameterise the signal 
traces, and the added dimension of detailed image timing information. The 
analysis described here is based on that which has been successfully employed 
for the analysis of Whipple 10 m telescope data for many years, with the use 
of image timing information being restricted to helping to define the optimum 
signal integration window. 

We have used a two-pass method in this analysis. In the first pass, a wide 
(20 ns) integration gate is applied to each FADC trace at a fixed time po- 
sition in order to calculate the integrated charge and the pulse arrival time. 
To (defined as the half-height point on the pulse leading edge). PMTs are se- 
lected which produce an integrated signal greater than 5 times their pedestal 
standard deviation, or greater than 2.5 times the pedestal standard deviation 
while also being adjacent to a PMT selected at the higher threshold. The re- 
sulting shower image is then parameterised with a second moment analysis, 
the results of which can be described by an ellipse as shown in Figure 10 [3]. 
Figure 11 shows Tq as a function of PMT position along the long axis of the 
ellipse for this event; the time gradient across the image is given by a straight 
line fit to this graph. In the second pass over the image, a shorter (10 ns) 
integration window is applied to each trace, but the location of the window is 
calculated using the expected position of the pulse according to the measured 
time gradient across the image. 

The standard candle of ground-based 7-ray astronomy is the Crab Nebula, 
and the sensitivity of any new instrument is best tested against this bench- 
mark. The first VERITAS telescope saw first light in February 2005, leaving 
only a small time window in which to collect data on the Crab while the 
source was still observable at high elevations. Observations were taken in an 
ON/OFF mode, wherein a 28 minute ON exposure at the source position is 
followed or preceded by a 28 minute OFF-source exposure, at the same el- 
evation but offset by ±30 minutes in right ascension. The total ON-source 
exposure was 3.9 hours with the same amount of OFF-source data collected 
in order to provide a measure of the cosmic ray background. Due to the purely 
electromagnetic nature of 7-ray initiated air showers, 7-ray Cherenkov images 
are expected to be more compact than cosmic ray images and can be prefer- 
entially selected by applying cuts based on the dimensions of the event ellipse 
{length, width) [3]. The parameter a describes the orientation of the long axis 
of the image ellipse relative to the line joining the centre of the field-of-view 



and the ellipse centroid. Showers originating from a point source, such as the 
Crab Nebula, will be oriented with their long axes pointing back to the po- 
sition of the source at the centre of the camera. The angular distance of the 
image from from the source position is related to the distance of the shower 
core position on the ground. Finally, the ratio of the image length to its total 
size (where size is the integrated charge over all PMT signals in the image) is 
used to discriminate 7-rays from the otherwise overwhelming background of 
local muons which generate short arcs in the camera with a constant size per 
unit arc length. Cuts on these various parameters were optimised using a data 
set of observations of the variable TeV blazar Markarian 421 when this source 
was known, from contemporaneous observations with the Whipple 10m tele- 
scope, to be in a high emission state. The resulting a histogram for the Crab 
Nebula observations is shown in Figure 12. The excess at low a corresponds 
to a significance of 19. 4o" (where a refers to one standard deviation above the 
background), indicating a sensitivity of ~ lOa for one hour of ON-source ob- 
servations. Figure 13 shows the significance map of the reconstructed source 
position using the method of Lessard et al. [23]. 

Following the Crab Nebula observations, the telescope was used in an ob- 
serving campaign throughout 2005, complementary to that of the Whipple 
10 m telescope, which resulted in > lOa detections of the known 7-ray sources 
Markarian 421 and Markarian 501 [24] and data sets on various potential TeV 
sources. A full online analysis package has been developed and was success- 
fully used to detect strong 7-ray flaring behaviour from Markarian 421 on a 
timescale of minutes. 



5 Simulations 



5.1 The Simulation Chain 



In order to accurately calculate source fluxes and energy spectra, it is necessary 
to develop a detailed model of the telescope performance, to which simulated 
air showers are presented such that the telescope detection efficiencies can be 
calculated. A complete chain of Monte Carlo simulations has been developed 
[25] . It consists of air shower simulations with CORSIKA [26] and a detailed 
simulation of the telescope response (GrISU [27]). CORSIKA version 6.20 is 
used with the hadronic interaction models QGSJet for primary energies above 
500 GeV and with FLUKA for primary energies below 500 GeV. Simulated 
hydrogen, helium, and oxygen nuclei-induced air showers were produced in 
an energy range from 10 GeV (150 GeV for oxygen) to 50 TeV for various 
telescope elevation angles. Spectral indices were taken from fits to balloon 
measurements [28]. The shower cores were chosen to be distributed randomly 



10 



within a circle of radius 1000 m centred on the telescope and the isotropic 
distribution of the cosmic ray incidence angles was simulated by randomizing 
the shower directions in a cone of radius 3.5° around the pointing direction 
of the telescope (while showers from outside of this cone angle can trigger 
the telescope, we estimate the fraction to be no more than 5% of the total). 
Primary 7-rays were simulated with a Crab-like spectrum {E""^'^) in an energy 
range from 10 GeV to 10 TeV on a circular area at the ground with radius 
450 m. 



Measurements of the atmospheric properties at the site of the first VERITAS 
telescope are currently in progress. The calculations described here used the 
U.S. standard atmosphere, which does not always reflect the properties of 
the atmosphere in Southern Arizona. Atmospheric extinction values were esti- 
mated using MODTRAN 4 [29] assuming 50 km visibility at 550 nm at ground 
level. The photoelectron rate per PMT is measured to be 100-200 MHz, cor- 
responding to a night sky background rate of ~ 2.8 x 10^^ photons m~^ s~^ 
sr~^, which was used in the simulations and produced a simulated pedestal 
standard deviation of a similar level to that observed in the data. 



The telescope simulations consist of two parts, the propagation of Cherenkov 
photons through the optical system and the response of the camera and elec- 
tronics. The geometrical properties of the optical system are fully implemented 
in the simulation, including allowance for the surface roughness of mirror facets 
and for random scatter in their alignment. The camera configuration corre- 
sponds to that of April 2005, i.e. a 499 pixel camera with 2.86 cm diameter 
phototubes without light cones, with a total field-of-view of 3.5°. The absolute 
gain calibration has been described in Section 3.1. In order to estimate the ex- 
act pulse shape to be used in the simulations, we examined pulses in the data 
produced by muon events and laser fiashes, both of which are light sources 
with a time profile which is expected to be shorter than the bandwidth of the 
electronics. The resulting average single photoelectron pulse has a rise time of 
3.3 ns and a width of 6.5 ns. Cherenkov photons hitting a PMT produce a sin- 
gle photoelectron pulse with appropriate amplitude and time jitters applied. 
Electronic noise and all efficiencies, including mirror refiectivities, geometri- 
cal, quantum, and collection efficiencies, and losses due to signal transmission 
were modeled. The pulses are digitized into 2 ns samples with a trace length 
of 24 samples refiecting the properties of the FADC system. The trigger sim- 
ulation uses a simplified model of the constant fraction discriminator, which 
is the first stage of the VERITAS multi-level trigger, and a full realization of 
the pattern trigger, requiring three adjacent pixels above threshold in a time 
window of 6 ns. The currently used trigger threshold of 70 mV corresponds to 
about 6.7 photoelectrons. The output of the telescope simulations, i.e. FADC 
traces for all PMTs, are written to disk in the VERITAS raw data format. 
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5. 2 Raw Trigger Rate 



The raw trigger rate of the telescope with trigger conditions as described 
above is ~ 150 Hz at high elevation. Dead time losses due to readout are 
~10%, the corrected trigger rate is consequentially ~160 Hz. The simulation 
of the cosmic ray background results in trigger rates of 101 Hz from air showers 
induced by protons, 26 Hz from helium nuclei, and ~5 Hz from nuclei of the 
CNO group for an elevation of 70°. The Monte Carlo calculations reproduce 
the observed trigger rate with an accuracy of ~20%, which is acceptable, 
taking into account an estimated uncertainty in the cosmic ray fluxes of ~25% 
and various systematic uncertainties in the modeling of the telescope and the 
atmosphere. 



5.3 Comparison with Real Data 



Simulations and real data were analysed in the same way. The analysis method 
used was a simpler, single-pass version of the analysis described in Section 4, 
which has a slightly reduced sensitivity. 7-ray candidates were extracted from 
observations of the Crab Nebula (3.9 hours ON-source) and Markarian 421 (4.2 
hours ON-source) in March and April 2005. All data were taken in ON/OFF 
mode at elevations above 60° and in good weather conditions. The 7-ray selec- 
tion cuts result in a combined significance of 21.7(T and 900 7-ray candidates 
for the combined observations of the Crab Nebula and Mrk 421. The measured 
7-ray rate from the observations of the Crab Nebula alone is 2.1±0.2 minute"^. 
Using the Crab spectrum reported in Hillas et al. [31] results in a simulated 
7-ray rate from the Crab Nebula of 2.2 minute"^. 

Ideally, the best image parameter cuts can be derived from the Monte Carlo 
simulations, which provide a data set of pure 7 rays with high statistics. 
Simulations also provide the only method by which to estimate the primary 
energy of the incident 7 rays. For the simulations to be useful, it is necessary 
to verify that the simulated images are an accurate representation of the data. 
Figure 14 shows the distributions of the image parameters a, width, length 
and distance for simulated 7 rays and for the 7-ray excess observed in the real 
data. The agreement is good, indicating that the simulations are valid. 



5.4 Crab Nebula Energy Spectrum 



The effective area of Telescope 1 after cuts is shown in Figure 15. The filled 
circles show the effective area for the hard 7-ray selection cuts described above 
which reject most of the high energy events due to the cut on distance. This 
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is a problem which is well known from our experience with the Whipple 10m 
telescope, and is solved by using size-dependent 7-ray cuts [30], which do not 
provide the same level of background rejection, but which take the dependence 
of the image parameters on primary energy into account and retain ~ 80% of 
the 7-ray signal. 

The maximum effective area after size-dependent cuts is ~ 2.8 x 10^ m^. The 
energy threshold of a Cherenkov telescope is conventionally defined as the 
position of the peak of the energy spectrum of the source convolved with the 
effective area curve of the detector. According to this definition, the threshold 
is 150 GeV at trigger level, 160 GeV after applying szze-dependent cuts, and 
370 GeV after applying the hard cuts described earlier. 

In order to reconstruct the energy spectrum of 7 rays from the Crab Nebula 
we have used the method outlined in Mohanty et al [30]. Size-dependent cuts 
are applied to the data, and the primary energy of the resulting 7-ray excess 
events is estimated using a function which depends on the image size and 
distance; the exact form of the function is derived from the Monte Carlo 
simulations. Figure 16 shows the resulting reconstructed energy spectrum of 
the Crab Nebula. A power-law fit to the data points gives a spectral index of 
2.6 ± 0.3 and a differential flux at 1 TeV of (3.26 ± 0.9) • 10"^ m-^s-^TeV"^ 
(statistical errors only). This agrees well with earlier measurements by other 
telescopes [31,32]. 



6 Discussion 



The first VERITAS telescope was operated throughout 2005, has met all tech- 
nical specifications and has detected a number of TeV 7-ray sources. The sen- 
sitivity, at ~ lOcr for one hour of ON-source observations of the Crab Nebula, 
is already appreciably better than the Whipple 10 m telescope (at ~ 7a for 
one hour of Crab observations); however, in order to provide the level of sen- 
sitivity required from this generation of detectors, it is necessary to reject the 
muon background at the hardware level. This can only be achieved by the in- 
stallation of further telescopes. The major mechanical components of all four 
telescopes have now been delivered to the Mt. Hopkins base camp. The second 
VERITAS telescope has recently been installed 85 m away from the first tele- 
scope on an East- West baseline. Figure 17 shows the length/ size histogram 
for the first telescope for all events which trigger, and for those which generate 
a trigger in both telescopes within 10 fis of each other. The distinctive peak 
due to local muon events in the histogram of all events is clearly removed when 
the two-telescope trigger is imposed. A true "hardware" two-telescope trigger 
is currently being installed, which will correct for changing time delays due 
to the source movement in the sky, and reduce the required coincidence time 
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to ~ 50 ns or less. Only events which satisfy the two-telescope trigger will 
be read out by the data acquisition system, reducing the data rate dramati- 
cally and so enabling us to lower the individual telescope trigger thresholds, 
and hence the energy threshold of the experiment. The angular resolution for 
point source analysis with a single telescope is 0.17° (Figure 13); additional 
telescopes will provide multiple views of each shower, reducing the angular 
resolution by around a factor of two. 

The agreement of the Monte Carlo simulations with observational data illus- 
trated here shows that we have a good understanding of the first VERITAS 
telescope and that the design performance is being met. This being the case, 
the predicted behaviour of the four-telescope array system is likely to be reli- 
able, indicating a flux sensitivity of 3.4 x 10~^^ cm~^ s~^ at 100 GeV for a 50 
hour exposure 
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Fig. 1. The VERITAS Telescope 1 as installed at the Whipple Observatory base 
camp. The collector dish has a diameter of 12 m and a focal length of 12 m and 
comprises 350 mirror facets. A 499-PMT camera is installed in the box at the 
focal point. The buildings in the foreground house the electronics and power supply 
systems. 
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Fig. 2. The azimuth and elevation residuals (difference between measured and 
requested position) for a short tracking run. 
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Fig. 3. A close-up view of the hexagonal VERITAS mirror facets. The three ad- 
justment points which are used in aligning each facet can be easily seen. 
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Fig. 4. Average mirror reflectivity as a function of wavelength for photons normal 
to the mirror. Measurements were made in the laboratory after the mirrors were 
produced. 
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Fig. 5. An image of Polaris in the focal plane of the telescope recorded with a CCD 
camera. The point spread function is 0.06° FWHM. The circle indicates the size of 
a VERITAS PMT (0.15° diameter). 




Fig. 6. The 499 PMT pixel camera. The focus box is 1.8 m square. A remotely 
operated shutter usually covers the camera during daylight hours. 
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Fig. 7. The trigger rate as a function of CFD threshold for two different topological 
trigger configurations. 
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Fig. 8. A FADC trace produced by Cherenkov light from a cosmic ray air shower in a 
single signal channel. The dashed horizontal line is the electronic pedestal level, the 
dashed vertical line shows the pulse arrival time, to, and the shaded area indicates 
a 10 ns (5 FADC samples) integration window. 
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Fig. 9. The single photoelectron response for one PMT at standard operating gain. 
The fit assumes a Poisson distribution of photoelectrons and a Gaussian distribution 
for the integrated charge produced by a single photoelectron. The average conversion 
factor over all PMTs is 0.19 ± 0.02 photoelectrons per digital count. 
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Fig. 10. The charge distribution across the camera for a cosmic ray event (the grey 
scale is in digital counts). 
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Fig. 11. The Cherenkov pulse arrival time distribution (in units of FADC samples 
= 2 ns) along the long axis for the cosmic ray image shown in Figure 10. A PMT 
position of zero corresponds to the centre of the image ellipse. Error bars correspond 
to the signal size-dependent time resolution for each PMT, as measured using the 
laser calibration system [17]. 
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Fig. 12. The image orientation angle, a, as described in the text, for ON and 
OFF-source observations of the Crab Nebula after 7-ray selection cuts. The excess 
events at low values of a are due to 7-rays originating from the direction of the 
Crab Nebula. Selecting events with a < 7° results in a detection with a significance 
of 19.40-. 
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Fig. 13. The two-dimensional map of reconstructed source position in degrees on 
the camera plane (the camera is 3.5° in diameter). The Crab Nebula was located at 
the centre of the camera. The angular resolution is 0.17° (68% containment radius). 
Note that adjacent bins are not statistically independent. 
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Fig. 14. Comparison of image parameter distributions from real data (closed sym- 
bols) and Monte Carlo calculations (open symbols). The parameters are described 
in the text. 
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Fig. 15. Effective area of Telescope 1 after hard 7-ray selection cuts (filled circles) 
and size-dependent cuts (filled triangles). 
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Fig. 16. Energy spectrum of the Crab Nebula compared with earUer results from 
the Whipple 10 m telescope [31] and HEGRA [32]. 
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Fig. 17. length/ size distributions (parameters defined in the text). The single 
telescope histogram includes all events which trigger the telescope. The distinctive 
peak in this distribution is due mainly to local muon events. The two telescope 
distribution includes only events which triggered both telescopes within a 10 /xs time 
window. Events due to local muons are removed when this condition is imposed. 
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